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1 Introduction 



We consider the stochastic boundary value problem of elliptic type 

Au(x) - f(u(x)) = g{x) + F(x), for x G D (1.1) 
u(x) = 0, for x £ dD, 

where D = (0, l) d , d > 4, and F( x) is the coloured noise defined at the 
beginning of Section |2 The coefficient g belongs to L 2 (D) and / is of the 
form 

f(x) = f 1 (x) + f 2 (x), (1.2) 
with /, fx, /2 : R — > M satisfying 

(/I) /i is continuous, non-decreasing and bounded, sup^g^ |/i(x)| < M, 
(/2) /2 is Lipschitz with small Lipschitz constant L. 

Let Ge> be the Green function of the Poisson equation Av = b, with 
boundary condition v (x) = if x £ dD. By a solution of (jl.l|) we mean a 
stochastic process (u(x),x S -D) satisfying 

u(x) = / G D (x,y)f(u(y))dy + / G D (x,y)g(y) dy + / G D (x,y)dF(y). 
Jd Jd Jd 

(1.3) 

For dimensions d = 1,2,3, -F a white noise and .D a bounded domain 
in R rf with regular boundary, the existence of a unique solution to equation 
(|1.3j) has been proved in [3] (see also |lj). In this framework, for the particu- 
lar domain D = (0, l) d , Gyongy and Martinez introduced in a numerical 
scheme based on discretization of the Laplacian and gave the rate of conver- 
gence in the L 2 (Z))-norm (see also 2 for related work). For these range of 
dimensions the Green function satisfies sup x&D \\Gd(x, -)\\l 2 (d) < °°- This 
fact has two important consequences. Firstly, (J D GD(x,y)dF(y),x G D) 
defines a Gaussian process; secondly, one can work with the Fourier series 
of Go- Actually, the lattice scheme in is related to truncated Fourier 
expansions. 

For parabolic equations driven by space-time white noise, numerical 
schemes based on lattice approximations have been introduced in jS] (see 
also ^01)- The method set up in these papers has been successfully applied 
in and also to other types of SPDEs, for instance, to a wave equation 
in spatial dimension 1 in |14j . 

In this paper, we study a lattice approximation for the stochastic partial 
differential equation (jl.lj) for d > 4, extending the results of |llj . 

For dimensions d > 3, the Green function associated with ()1.1|) is given 

by 

G D (x,y) = G(x,y) + E x (G(B T ,y)), with G(x,y) = _° d (1.4) 
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In these expresions, Cd is a constant and B is a d-dimensional Brownian 
motion starting at x, stopped at r -its first exit time of D (see |S] for 
details) . 

We will prove that for d > 4, G D (x, •) G L a (D), a G [1, d/(d - 2)), uni- 
formly in x. Hence, we cannot use a L 2 theory. In particular, the stochastic 
integral J D GD(x,y)dF(y) with respect to a white noise cannot be defined 
as a real-valued I? random variable. This problem forces the choice of a 
coloured noise, as a way to give a rigurous meaning to Ijl.ljl . 

The contents of the paper are as follows. Section [21 is devoted to study 
the spde. First, we attach a precise meaning to the stochastic integral term 
in Q1.3JI as a Gaussian process indexed by elements of its reproducing kernel 
Hilbert space and give sample path properties related with the regularity of 
the covariance of the noise. Secondly, we prove existence and uniqueness of 
solution. The methods are common to those of nonlinear monotone operator 
equations (see ^2], ^1) an d those used also in [3] and @]. However, since 
our setting is different, we feel interesting to give the details. We also ana- 
lyze properties of the solution; in particular, Holder continuity of the sample 
paths. Section |3] is devoted to the lattice approximation and the analysis 
of the rate of convergence (see Theorem 0J). First, we introduce a family 
of smoothed Green functions obtained by convolution of Gd(x, •) with an 
approximation of the identity. Then, we introduce a lattice approximation 
which corresponds to a weighted discretization of the Laplacian on a given 
grid, with weights related to the smoothing. The approximation result fol- 
lows from a careful analysis of the errors in both of these approximations, 
taking as smoothing parameter an appropriate function of the norm of the 
grid. 

Let us give some indications about notation. By \x\ we denote the 
Hilbert-Schmidt norm of any element x G M. d . The letter C denotes, un- 
less otherwise stated, a constant that may not be the same from one oc- 
currence to another. Sometimes, we denote fixed values of constants by 
adding a subindex, and the dependence on some parameters (as the di- 
mension) with an argument, e.g. C\ or C(d). Throughout the work, the 
symbols (3 = (/?i, . . . ,/3rf), i = («i, ■ ■ ■ ,id) denote indexes belonging to the 
sets I d = {1, 2, . . . } d , I d = {1, . . . ,n-l} d , with |/?| 2 = /3 2 + - • -+/3j. Observe 
that for (3 G I d \I d , \/3\ 2 > n 2 , and for (3 G I d , d < \(3\ 2 < dn 2 . 

2 Study of the equation 

Let if be the density of a non negative measure on W d , non negative definite 
and tempered. We consider a centered Gaussian process F(ip), indexed by 
the space T>(M. d ) of Schwartz test functions, with covariance function 




(2.1) 
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ifo.,ih€ V(R d ), defined on some probability space ($7, J 7 , P). 

Let Ti denote the completion of the inner product space consisting of 
functions ip £ T>(R d ) endowed with the inner product 

(iI>i,^2)h= / tpi(x)(p(x-y)tl) 2 (y)dxdy. (2.2) 

JR d JR d 

The space Ti is the reproducing kernel Hilbert space corresponding to F. 

We consider the stochastic partial differential equation (jl.lj) . with the 
assumptions given in Section ^ We notice for further use the following 
property: 

(P) / is a function of the form (jl.2j) with f\ non-decreasing and Lipschitz 
with Lipschitz constant L, if and only if for every u, v £ R, 

(u-v)(f(u)-f(v))>-L(u-v) 2 . (2.3) 



2.1 The stochastic integral 



In order to give a rigourous meaning to (jl.lj) . we have to precise what type 
of stochastic integral we are considering. By classical results on abstract 
Wiener spaces, h £ Ti i— > F(h) defines a linear continuous functional that 
satisfies (|2.1j) (see HUD- Thus, the stochastic integral term in (|1.3j) is well 
defined as long as we prove that Gd(x,-) £ Ti. The next lemma provides a 
sufficient condition for this property to hold. 

Lemma 1 Let p £ [1, oo), p + ^ = 1. 

1. (1) Let ip £ L p (R d ), ip £ L^(R d ). Then, 



ip(x)tp(x — y)ip(y) dx dy 



< 



2. (2) Assume V? £ L 1 (R d ) n 



),<f £ LP 



,4>{x)y(x - y)ip(y) dx dy 



< 



LP 



LP(U. d )W\\ ^ 



Then, 

L 1 {R^W^W Lp' 



Proof. Applying first Holder's inequality and then Young's inequality for 
convolutions (pQ, Corollary 2.25) yield 



ip(x)tp(x — y)ip(y) dx dy 



■ip(x)(ip * ip)(x) dx 



< 



LP (S. d )W 



\LP' 



< 



LP(M. d )W ( p\\LP'/ 2 



Young's theorem (pP, Theorem 2.24) implies 



ijj(x)(ip * ip)(x) dx 



< 



L 1 (R d )W\\LP' 



LP 



Both estimates yield the lemma. 

We next prove a basic result on the Green function Gd- 



□ 
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Lemma 2 For any a G [1, ^2), i/iere exists a positive constant C\ depend- 
ing on a and d, such that 

sup \\G D (x,')\\ La rm < Ci. 

x&D 



Consequently, ||Gd||l<*(Dx.D) < G. 
Proof. Set h = \\G(x, -)||z«(.D)- Clearly. 



, , < C{d) / r d(l-«)-l+2« d r . 



The last integral is finite if and only if a < ^2', in this case, its value is a 
constant C(a,d) independent of x. 

Let I2 = \\E X (G(B T , y))\\L a (D) an d let P x be the law of the random 
variable B T . By Minkowski inequality, 



G(u, •) dP x (u) 



< 



L a (D) 



\\G(u,-)\\ La(D) dP x (u)<C(a,d). 



Since bv (|1.4j) . ||Gd(x, -)l|.z/*(.D) < /1+/2, with the upper bounds obtained 
so far, we finish the proof of the lemma. 

□ 

Fix a G [1, 3^2) and denote by a' its conjugate. Set C a = L a ' (R d ) U 
L a / 2 (M d ). The preceding lemmas show that if (p G £ a , then, for any a; G R rf , 
Gd(x, •) G H. Consequently, the stochastic integral J^, Gd(x, y)dF(y) is well 
defined and 



E 



r 2 
j D G D (x,y)dF(y) < \\G D (x r )f La[D) < C\. (2.4) 



In addition, for each p G [l,a'], Holder's inequality along with the hyper- 
contractivity property and (j2.4|) yield 



E 



D 



G D (x,y)dF(y) 



L a '(D) 



1/p 



< E 



D 



G D (x,y)dF(y) 



L a '{D) 



1/a' 



< ( '[ J I I. 



D 



G D (x,y)dF(y) 



2 X a'/2 N 1/a' 

dx 



< 



C^J D \\G D {x,-)\\i a{D) dx^ ' <CC 1 



(2.5) 



We next prove that the stochastic integral defines a Holder continuous ran- 
dom field. 
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Theorem 1 Fix A G (0,1). Assume that tp £ C a (D) with a in the inter- 
val (1, ((j_2)(2-A)v(d-i)A )• Then, for any d > 4, i/ie Gaussian random field 
{v(x) = J D GD(x,y) dF(y),x £ D} satisfies 

E(\v(x) -v{z)\ 2 ) < C\x-z\ 2X . (2.6) 

Therefore, a.s. the sample paths are Holder continuous of order 7 £ (0, A). 

PROOF. Fix x,z £ D. By the first inequality in (|2.4|) . 



E 



D 



G D (x,y)dF(y)- / y) dF(y) 



D 



<T 1 +T 2 , 



with 



T 1 = \\G(x,-)-G(z,-)\\l a{D) , 

T 2 = \\(E x (G(B T ,y))-E z (G(B T ,y))\\l a{D) . 

Schwarz inequality implies 



T 2 



\\G{x,-)-G{z,-)\\l HD) <C 
1 1 



D 



\x - y\ 



D 



\x - y\ 



d-2 



\ z - y\ 



d-2 



a(2-A) 



d-2 

1 
\ 2 



\ z - y\ 



d-2 



a A 



dy 



dy 



The last integral is finite provided that a & (0, tt "' 

For the first one, we apply the mean value theorem to obtain 



(d-2) (2- A) - 



D 



1 



1 



|x — y\ d 2 \z — y\ 



d-2 



oA 



dy < \x — z 



aX 



D 



d-3 



|x — y| d 2 |z — j/ 1 



d-2 



aX 



dy, 



with £ = fj,(x — y) + (1 — /x)(z — y) for some \x G (0, 1). Thus, if we define 
Di =Dn{\x-y\ < \z-y\} c 



< 



dy 



| x _ y^d 2| z _ y|d 2 

2|z-y| d " 3 " A 



< c 
c 



x — y\ d 2 \z — y\ d 2 
dy 

Di k - yl^" 1 ) 



+ 



D 2 



dy 



< C 



| x _ y\aX(d-l) — J Q r aX{d-l) 



= D\Dt 


we 


have 


dy+ [ 

Jd 2 

dy 




2\x - y\ d ~ 3 


\x - 


- y\ d ~ 2 \z — y\ d ~ 2 

) 


\ z _ y\aX(d-l) 
r d-l ^ r 
j.a\(d—l) 





aX 



dy 
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This last integral is finite if and only if a G (0, jjzjjx)- Summarising, for 

2)(2-A)V(d-l)A h 

\\G(x,-)-G(z,-)\\ La{D) <C\x-z\\ (2.7) 



a ^ C 1 ' (d-2)(2-AW(d-lu) 



For the second term, let us observe that by the strong Markov property, 
E z (G(B T ,y)) = E X (G(B T — x + z,y)), and therefore 



(D) 



T 2 = C\\E x (G(B T ,y))-E x (G(B T -x + z,y))\\ 2 La 

< c(jj(G{u,-))-G{u-x + z,-)\\ La{D) dP x {u)\ <C\x-z\ 2X , 

where in the last inequality we have applied 1)2.7)) . Thus, 

\\G D (x, •) - G D (z, -)\\ LHD) < C\x - z\ x . (2.8) 

Since the process (v(x),x G D) is Gaussian, the statement follows from 
Kolmogorov's continuity criterion. 

□ 

Remark: Assume p G C a for some a G (1, jtt)- A slight variant of 
the proof of Theorem ^ yields Lipschitz continuity of the sample paths of 
(v(x),x G -D), a.s. 



2.2 Existence and uniqueness of the solution 

We give in this section a theorem on existence and uniqueness of solution 
for equation and we also establish some of the properties needed later. 

We shall often use the next property proved in Lemma 2.4; it is a 
consequence of the solvability of the Dirichlet problem on D and Poincare's 
inequality (see jH]). 

(P) There exists a constant a > such that for any ip G L 2 (D), 

Id ( y f D GD ^ x,y " >ip ^ dy ) v(. x ) dx - ~ a J D ( K J D G D{ x ,y)^(y)dy^ dx. 

(2.9) 

Unless otherwise stated, along this section A is a fixed number in (0,1), 
a G [1, (rf_2)(2-A)v(d-i)A )' a ' ^ s ^ s conjugate, that is, ^ + ^ = 1 and a the 
positive constant given in ()2.9)) . 

Theorem 2 Suppose that g G L a ' (D), f is of the form U.fy) and satisfies 
(fl) and (f2) with a Lipschitz constant L < min(a, C-f ). Assume also 
that the correlation density tp belongs to C a . Then, there exists a unique 
stochastic process solution to 
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Proof. Let 

B = {w : we C(D), w\ dD = 0}, 
and consider the operator T : B — > £>, denned by 

T(w)(x) = w(x) - / G D (x,y)f(w(y))dy. 



D 



Clearly, B C L a '(D). 

By Holder inequality and (|2.8|) . 



(G D (x,y) - G D (z,y))g(y) 



^ l|5lli«' (0)11^(37,1/) - ^-dO^Hl^D) 

< C|x-z| A . 
Together with Theorem^ this implies 

6(x)= / G D (x,y)g(y)dy+ f G D (x,y) dF(y) G B, 

JD JD 

for each u, a.s. 

We next show that the operator equation Tw = b has a unique solution 
for any b 6 B, which happens if T is a bijective operator on B. Uniqueness 
guarantees the measurability of the process (w(x),x £ D). 

Let us first check that T in one to one. Fix u and v such that Tu = Tv. 
Then, u(x) - v(x) = J D G D (x,y)(f(u(y)) - f(v(y)))dy. Multiplying both 
sides of this equation by f(u(x)) — f(v{x)) integrating over D and using 
we obtain 



D 



(u(x) - v(x))(f(u(x)) - f(v(x))) dx 



D (f(u(x))-f(v(x))) ^J D G D (x,y)(f(u(y)) - f(v(y)))dy) dx 

2 



< -a f (f G D (x,y)(f(u(y)) - f{v{y)))dy\ dx 
= —a I (u(x) — v (x)) 2 dx. 

D 



By 



(u(x) — v(x))(f(u(x)) — f(v(x))) dx > —L / (u(x) — v(x)) dx. 

D JD 

Hence, (a — L) f D (u(x) — v(x)) 2 dx < 0. Since L < a&ndu,v G L a ' (D) with 
a' > 2, this implies u{x) = v(x) for almost every x £ D. 

We next prove that T is onto, proceeding in a similar way than in [3]. 
In the next arguments, w is a fixed element on a set of probability one. 



7 



Step 1. A solution for a regular problem. Let b G £>, and b n G C^°(D), n > 1, 
such that 6 n — > 6 in L Q (£>). Clearly, the convergence also holds in L 2 (D). 
We construct in the next LemmaEla sequence of functions solving Tu n = b n 
such that u n — > it in L 2 (D); the limit -u will be our candidate for solution. 

Let us recall a basic result on the solution of nonlinear monotone operator 
equations (see e.g. |12l Theorem 2.1, pg. 171], or alternatively 16., Theorem 
26.A, pg. 557]): 

(E) Let X be a reflexive Banach space; denote by X* its topological dual. 
Let B : X — > X* be a strictly monotone, coercive, hemicontinuous operator. 
Then, for any G X*, the equation Bw = k has a unique solution on X. 



Lemma 3 For every n > 1, t/ie boundary value problem 

Au ra (x) — f(u n (x)) = Ab n (x) for x £ D, u n (x) = for x G 3-D, 



/ias a unique solution u n (x) G Wq' 2 (-D). 



Proof of Lemma EI Set X = W 1,2 (-D), and define 5 : Wo' 2 (£>) 



4,2/ 



(W 1,2 (£>))* by £u = -A« + /(«); that is, for any w G Wo'^-D), 
(Bu,w) = / Vit;(x) • Vu(x) + / w{y)f{u(y))dy 



D 



D 



The assumptions on / imply that this operator satisfies the properties re- 
quired in (E). Moreover, k := -Ab n G Cf{D) C (Wq' 2 (D))*, for any n > 1. 
Thus, the lemma follows from (E). □ 

The sequence {u n ,n > 1} satisfies 

u n (x) = G D (x,y)f(u n (y))dy + b n , for x G I?, u n | 9jD = 0. (2.10) 

JD 

Let us check that it is a Cauchy sequence in L 2 (D). By multiplying both 
sides of the equation satisfied by u n (x) — u m (x) by f(u n (x)) — f(u m (x)), 
integrating over D, and owing to (|2~3|) and (|2~9|) . we obtain 



Il2(d) 



+ a 



D 



2 

L 2 (D) 



< 



D 



G D (x, y)(f(u n (y)) - f{u m {y))) dy 
(b n (x) - b m (x))(f{u n {x)) - f(u m {x))) dx. 
By using ()2.1(jj) . and the fact that each u n G L 2 (D), we obtain 

(a — L)\\ll n — UmW^tD) + a ll^ri ~~ ^m|li 2 (D) 

< / (b n (x) -b m {x))[{f(u n (x)) - f(u m {x))) + 2a(u n (x) - u m (x))]dx 

JD 

< \\K - b m \\L 2 (D)[ 2M + L\\ u n ~ u m \\L 2 (D) + 2a||n n - u m \\ L 2( D )}, 
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that is, 

\\u n — u m \\ 2 L 2( D ^ < C\\b n — b m \\ L 2^ D -j[l + \\u n — u m \\ L 2( D )}, 

which implies that \\u n — w m ||L2(D) — * as n,m -> oo. 
Set u = \\m. n u n in L 2 (D). 

Step 2. u is the solution. We must prove that u G B and verifies (|1.3[) . That 
is, we would like to take limits in (12.10(1 . 

We choose subsequences u n and b n (still denoted with the same sub- 
scripts) converging to u and b almost everywhere. We proceed in three 
steps. 

Step 2.1. Assume that / is bounded (and continuous). Then 

u[x) - I G D (x, y)f(u(y)) dy = b[x) for x G D, u\ dD = 0, 

JD 

by bounded convergence, and u G £>. 

Step 2.2. Assume that / is bounded from below, that is, f(x) > — iV for 
every x and some > 0. Set f n (x) = fi(x) + (/2(a?) A n), n > 0. Notice 
that each / n satisfies (/l) and (/2). Let 

u n (x) = / G D (x,y)f n (u n (y))dy + b(x) for x G D, n n | 9D = 0, 
JD 

be the solution constructed in Step 2.1. 

We will need the following comparison statement. Its proof is very similar 
to that of Lemma 2.6 in and therefore omitted. 

Lemma 4 Let f and h satisfy (/l), (/2) and f(x) > h(x) for every x G R. 
Lei 6 G L Q ' (D) and 

u(x) - \ G D (x,y)f(u(y))dy = 6, 
w(x) - / G D (x,y)h(v(y))dy = b. 

JD 

Then, u(x) < u(x) /or almost every x £ D. 

The sequence (f n ,n > 0) is increasing; hence, by LemmaEJ the sequence 
of functions (u n , n > 0) satisfying 

u n (x) - G D (x,y)fi(u n (y))dy - G D (x,y)(f 2 A n)(u n (y)) dy = b(x), 

JD JD 

(2.11) 

is decreasing. Set u(x) = mi n u n (x). Notice that it is an a.s. finite function. 
The function fx being bounded, we can take the limit inside the first integral 
in the left hand-side of (|2.11j) . It remains to prove that we can also take the 
limit inside the second integral. For this, we need some a priori estimates 
provided by the next two statements. 
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Lemma 5 The sequence (u n ,n > 0) defined in \2.11\) satisfies 

SUp || U n || La i < ft, 

n v ' 

(A/+|/ 2 (0)|)C7 1 + |16i| La , 

ioze/i k = i-LC a ^i = su Pxe-D II^M^ 'JllL a (-D) ( see Lemma 

n. 

Proof. Since / is bounded from below, f n is bounded by some constant 
depending on n and thus each u n E L Q (D). Now, K/2 A n)(it n (a;))| < 
|/2(0)| + L\u n {x)\. Thus, by Holder inequality 

\u n (x)\ < (M + |/ 2 (0)|)d + LCiKH^pj + l&(*)l, 
and by integration on D, 

\\ u n\\ L a ' (D) 

< (M + |/ 2 (0)|)Ci + LCiKH^^) + \\Hl«'(d)- 
This yields the lemma. □ 
Corollary 1 | Ml < K, with k as in Lemma\^ 
Proof. It is a consequence of Fatou's lemma. In fact, 

I \u\ I £,«'(£>) < hminf ||ii n || LQ /^ D) < sup ||u n || LQ /( D) < k. 

□ 

Since u n {x) [ u(x) > —00 almost everywhere, we have |^(^)IX{ u >o} — 

\ U o( x )\X{u>0}i aIld \ U ( X )\X{u<0} ^ su P n K(x)|X {tt<0} . ThuS ' 

I tin (a?) I = l«n(aOIX{u>0} + \ u n(x)\X{u<0} 

< \Mx)\X{u>0} + \ U ( X )\X{u<0} = <P( X )> 

with (p E L a> (D) and therefore, 

|(/ a An)M*))| < l/aM*))l <M+ |/2(0)| +L^(i) E L Q '(L>). 
Thus, by the dominated convergence theorem and the continuity of f 2 , 



lim / G D (x,y)(f 2 An)(u n (y))dy = G D (x,y)f 2 (u(y))dy. 

n -* 00 JD JD 

Summarizing, if / is bounded from below, there exists u satisfying 
u(x)= G D {x,y)f{u{y))dy + b(x) for x E D, u\ dD = 0, 

JD 



and [|tt[| < k. 
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Step 2.3. / satisfies (fl) and (f2). Set f n = fi + (/2 V (-n)), n > 0. By the 
results obtained in the previous step, there exists u n such that 

u n {x) = / G D (x,y)f n (u n (y))dy + b(x) for x G £>, u| 9D = 0, 



and sup n ||n n || LQ / (D) < k. 

The sequence (/ n ,ra > 0) is decreasing; hence, by Lemma|l{u n ,n > 0} 
is increasing. Set u(x) = sup n « n (a;) for a.e. x. As in Step 2.2, it suffices to 
prove that 

lim / G D (x,y)(f 2 V (-n))(u n (y)) dy = G D (x,y)f 2 (u(y))dy. 

n ^°°JD JD 

Lemma El and Fatou's Lemma implies || HI £<*'(£)) < k. 

Since u n {x) | u(x), we have \u(x)\X{ u>0 y = sup n \u n (x)\X{ u>0 y, and 
H x )\X{u<0} ^ \Mx)\X{u<o}- ThuS > 

\u n (x)\ = \u n (x)\X{u>0} + \ U n(x)\X{u<0} 

< H X )\X{ U>0 } + \M x )\X{ u <o} = 

with V € L a '(D). 

Observe that \f 2 V (— n)| < I/2I. Therefore 

|(/ 2 V (-n))K(y))| < |/ 2 («„(y))| < M + |/2(0)| + LV(*)- 
Thus, by the bounded convergence theorem and the continuity of f 2 , 



n— >oo 



lim / G D (x,y)(f 2 \/ (-n))(u n (y))dy = G D (x,y)f 2 {u(y))dy. 

~ 'D ./D 



Hence, we have proved the existence of u G L a ' (D) satisfying 

1)= / Gd(x, y)f(u(y)) dy + 6(x) for x £ D, (2.12) 



ii|<9L> = 0. The terms in the right hand-side of (|2.12|) belong to B; therefore, 
so does u. 

□ 

2.3 Properties of the solution 

The solution of (|l.lj) given in Theorem [21 possesses some important proper- 
ties, as follows. 

Lemma 6 With the same assumptions as in Theorem^ for any p £ [l,a'], 

IMIz>(fl;i/*'(£>)) - C- 
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Proof. The solution u to equation satisfies HitH^a'm) < K, where k is 

as in Lemma|SJ Hence, it suffices to check that b S LP (O; L a ' (D)) , for any 
p € [1, a'], where 



b(x) 



GD(x,y)g(y)dy +/ G D (x,y)dF(y). 
d Jd 



The function x ^ J D Gd(x, y)g(y) dy is continuous and deterministic; conse- 
quently it belongs to L p (0; L a (£))). By virtue of ((231) the same conclusion 
holds for the stochastic integral J D GD(-,y)dF(y). □ 

Theorem 3 IFii/i i/ie same hypotheses as in Theorem^ for any p £ [l,a r ], 
the solution u to equation (jl.lj) satisfies 



E (\u(x) — u(z)\) p < C\x-z\ pX . 



(2.13) 



Consequently, a.s. the sample paths are j-Holder continuous with 7 G (0, A- 
A) 



Proof. Clearly, (E(\u(x) - u(z) 



,i/p 



<Ya=iI{x,z), with 



h(x,z) 
h(x,z) 
h{x,z) 



E 



E 



E 



D 



(G D (x, y) - G D (z, y))f{u{y)) dy 



P\ i/p 



D 



{G D (x, y) - G D (z, y))g{y) dy 



D 



(G D (x,y)-G D (z,y))dF(y) 



By Holder inequality, Lemma El and (|2.8|) . we have 

h(x, z) < \\G D (x, •) - G D (z, -)IU«(d)(M + |/ 2 (0)| + LE(\\u\\ P La , {D) ) 1 ^) 
< C\x-z\ x . 

A similar bound holds for the term l2(x, z). 

The hypercontractivity property and (|2.6|) yield I^(x,z) < C\x — z\ x . 
Thus, we have proved ()2.13|) . The statement about Holder continuity follows 
from Kolmogorov's criterion. □ 



3 Numerical approximations in L 2 (D)— norm 

This section is devoted to give a finite differences scheme for the spde (|1.1|) , 
We have shown that in dimension d > 4, the Green function does not belong 
to L 2 (D). Thus, the method used in for d < 3 does not apply directly. 
Instead, we combine a smoothing of Gr> with a weighted discretization. 
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3.1 Construction of the smoothing G e D 

We start by introducing some technical background (see 0)- Let T d be 
the d-dimensional torus endowed with the Haar measure dx. Any function 
/ : [-1, l) d — > R is identified with F : T d — ► R defined by 

F(e inx ) :=F(j™ 1 ,...,e i ™ d ), 

that is, F(e i7TX ) = f{x). Then 




Consider the odd extension of Gd(x,u) m the y-variables, that is, if yj G 
(0, 1), j = 1, . . . , d, we define 

Gd(ui, - ■ ■ ,-yi, ■■■■,yd) = -G D (yi, ■ ■ ■ ,yi, ■ ■ ■ ,yd)- 
The new function, still denoted by Gd(x, •), is now defined on (—1, l) d . Let 

G x D (e^) = G D (x,y) 
its identification in the torus setting. Observe that Gp satisfies 

/ \G x D (e^)\ a dy = [ \G D (x,y)\ a dy = 2 d \\G d (x,-)\\ a La{D) . 

Let ip(x) G C£°((— 1, 1)) be an even function, < ip < 1, J_ x ip = 1 and 
= UliV^)- Clearly, G C™((-l,l) d ) and it is an even function 
in each variable X{. Set 

d d 
$(e*™) : = U^') :=I]^) = *(x). 

i=l i=l 

The functions 

e > 0, define an approximation of the identity in T d . 

By means of <3? £ (e l7r '), we define the smoothing of Gd(x, •) as follows: 

G% £ (e^) = / G x D (e i7r{y - u) )$e(e inu )du. 
Ji d 

Clearly, Gf^e^) G L a (T), therefore G x £ -► G X D in L a (T) as e -> 0. Define 

Gf ? (x,y) = G^(e^). 
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It is not difficult to check that, for any y = (yi, . . . , — y,, . . . , yd), i = 1, . . . ,d, 
we have G £ D (x,y) = —G £ D (x,y), that is, G £ D (x,y) is odd in all the variables 
yi. In addition, for any a 6 [1, ), 



sup sup \\G £ D (x,-)\\ La{D) <CCi, 

e x£D 

with the same constant C\ defined in Lemma El 

The next result provides a bound for the error of the smoothing. 



(3.1) 



Lemma 7 Fix A £ (0,1) and as in Theorem^ assume that ip £ C a (D) 
with a & (0, (rf_2)(2-A)v(rf-i)A )• There exists a constant C such that, for 
every e > 0, 

\\Gd - G e d\\l°'(dxD) < Ce x . 
Proof. Since Gd(%, •), G £ D (x, •) are odd in the y-variables, we have 
\\G d (x,-)-GUx,-)\\l«(d) = 2- d \\G D (x,.)-GUx,-)\\L«((-i,iy) 



d \\Gf) - & X jf\\ L <x(jdy 



Thus, by Fubini's theorem 

W^D - ^D\\L a {DxD) — Z I 



L a (D) llL"(T d )' 



with the variable x integrated over D. Holder's inequality with respect to 
the finite measure on T d given by <& £ {u)du implies, for any fixed y, 



j*( e «w)_G£ B ( e «Tr\ll« 



l x D (e l7Ty ) 



L a (D) 

Ue in(y - u) )Mu)du 



L a (D) 
a 



(G x D (e i7Ty ) - G x D {e l7riy ~ u) ))$e{u) du 



L<*{D) 



< 



lUe i7Ty )-GUe i7T{y - u) )f ra (n Mu)du. 



By (|2.8|) and Fubini's theorem, we see that 

\\Gd - Gh\\ a La(DxD) < C \u\ aX <S> e (u)du. 

Taking into account that & £ { u ) 7^ if and only if u G (— £,e) d , we conclude 
the proof. □ 
Next, we give the Fourier expansion of G £ D (x,y). Let us recall that for 
any d > 1, the set of functions 

vp(x) = sin(/?i7rxi) • • • sm.((3 d itx d ), [3 £ I d , 

is an orthogonal complete system in L 2 (D), with Ht^H^m) = 2~ d / 2 . 
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Lemma 8 For any e > 0, we have 

G e D (x,y) = ~¥T^~ vp{x)vp{y), 

in L 2 (D x D) and a.e., where denotes the Fourier transform of ^ . 

In the sequel, we shall often use the following remark: The function 
is rapidly decreasing, that is, for any multiindex 9 G (0, oo) , there exists a 
constant C(9) such that |£Htf(£)| < C(6). 

Proof of Lemma [HI By virtue of Young's inequality for convolutions 

\\GUxr)\\LHD) = 2- d \\G^(e i n\\ LW 

<2 ||Gfj|| ia ( T d)||$ e ||£ P ( T d) (3.2) 

= \\GD(x,-)\\ La{D) \\^e\\ L P((-l,l)d), (3-3) 

with \ + 1 = i + i. For any p > 1, sup e ||*e||ip((_i,i)<«) < C- Hence, by 
Lemma |21 

sup sup \\G £ D {x,-)\\ L 2 iD) < Ci, (3.4) 

and consequently Gf} E L 2 (D x D). 

We next compute the Fourier coefficients of Gfj, as follows. Set vp(y) := 
ap(e my ). By Fubini's theorem 

G £ D (x,y)vp(y)dy = 2- d / G e D (x,y)vp(y) dy 

- d [ a, x j £ (>y\ 



2- a / GS £ (e^)^(e^)dy 

T d [ ([ G^ie^y-^^eie^du^ade^dy 
J T d \J T d J 

T d f (f GUe iAv ~ u) )o-3(e i7ry )dy)$ £ (e i7ru )du, 
J T d \J T d J 



Set K(e i7TU ) = f Jd Gf ) (e i7r ( y - u ))cj / 3(e^) dy. By a change of variables, 

K(e i7TU )= [ G x D {e inz )o-a{e i ^ z+u) )dz= [ G D (x,z)v (z + u) dz. 
JT d J(-i,i) d 

We compute this last integral using the formula 

d 

vp(z + u) = J|sin/3j7r(zi + m) 

i=l 
d 

= J^[(sin/3j7T2;j cos/3j7rnj + cos/3j7TZj sin/3j7rtii). 
i=i 
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The product on the right-hand-side of the former inequality consists of a 
sum of terms, each of them being a product of d factors either of the form 
sin Piirzi cos fiiirui or cos Piirzi sin/3j7TUj. Since Go(x,z) is odd in all the 
z- variables, integrals of terms containing factors cos a-iirzi sin a^Ui will be 
zero. Thus, 

d 



K(e tnu ) = / Gd(x,z) T\ (sin PiTTZi cos j3iirui)dz 
IT cos /3iirui ) / G D (x,z)va(z)dz 



cos PiTTu^j 2 d J Gd(x, z)vp(z) dz 

J_|cos/3i7rnj )—2^2 v p( x ), 

where in the last equality we have used the properties of Go- Therefore, 
[ GUx,y)v (3 (y)dy = 2- d [ K(e^ u )^ £ (e m ) du, 

JD JT d 

= 2 d I ( ]J COS Pi^ u i]-2^H V p( X )^e(u)du 

^(-u) J v 1=1 ) 

d 



i\1\2 v p( x ) ft / cos /3iTTx^ £ (x) dx, 



But 



[ cos f3iirxil) £ (x) dx = [ e t/3t7VX dx = ip(eoci), 

J (-1,1) J(-e,e) e e 

where ip stands for the Fourier transform of ip in M. Therefore, 

f — 1 d - —1 

J^G £ D (x,y)vp(y) dy = -^^ V p{x) J[ V»(e/?i) = -^^vp{x)V{ep). 

This finishes the proof of the lemma. □ 

3.2 Construction of the discrete approximations 

Let us introduce some notation. For any n > 1, we consider the grid of 
D = [0, l] d given by 

£={- = (-,...,-): j fc = 0,l...,n, fc = l,..., d} CD. 
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For any point ^ £ Q, we set = 

n n (x) = ^ for each x £ Dj. 

On the space X = {u : u = {ui} ieI d} = R( n_1 ) endowed with the 
Hilbert-Schmidt norm, we consider the second order difference operator A : 
X -> X defined by 

d 

(Au)i = ^2n 2 [u^ ej - 2m + u i+ej ] , 

3=1 

where {ej} d =1 is the canonical basis of M. d , The set of vectors of X, 

^) d/2 U,^Glf}j,(U,) t =v^), l Gl d n 

is an orthonormal system in X of eigenvectors of A, with eigenvalues 
A^ = -tt 2 (/3 1 % 1 + ---+^ % ), 

where q = sin 2 (j^J (i^) • Notice that 4j- < q < 1 for every 1 < I < 
n — 1. 

In the sequel, we consider the lexicographic order in I d . Denoting by U 
the (n— l) d matrix whose rows are the vectors Up j: (here j3j, j = 1, • • • , (n — 
l) d , denotes the lexicographic enumeration of I d ) we have 

A = U l DU, 

with D the square diagonal matrix with entries Djj = X^. 

For any e > 0, define D e the square diagonal matrix in dimension (n — l) d 
with diagonal elements 

We also consider a sequence (g n ,n > 1) of step functions defined on D 
such that for n big enough, — <7 n || 2 a ,^ < C/n 2 for some positive constant 

C. 

The discrete approximations of u are defined first on points of Q as 
follows. If ^ £ PI <9-D, define u^(^) = (boundary conditions). For ^, 
with j £ I d , we define "U^(n) to be the solution of the system 

(U'D'U)^ = /«) + 5n + n d F, (3.5) 

where F is the vector (F(Di),i G Z^), with = X D .d-F(y), and 

5n(o;) = g n (i<-n(x)), n > 1. Finally, for any x £ D we define = 
<(k„(x)). 



11 

n 



,h+i 

n 



n 



and define 
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We prove in Proposition^that a solution to equation Q3.5J) exists. More- 
over, proceeding as in it is easy to check that u £ n satisfies the mild 
equation 

u £ n (x) = G £ D (x,y)f(u £ n (y))dy + G £ D (x,y)g n (y) dy (3.6) 
Jd Jd 

G £ D , n (x,y)dF(y), (3.7) 

D 
with 

G D , n (x,y) = \ v/3(K n (x))vp(K n (y)). (3.8) 

fie* Xf> 

We will prove later that an appropriate sequence u n := Un"^ of these ap- 
proximations converges to the solution of (jl.lj) in the space L p (O; 1? (D)) , 
for any p £ [1, a'}, with an specific rate. 

Proposition 1 Assume that ip £ C a for some a £ [1, ^4). Suppose also 
that g £ L a ' (D), f is of the form il.ty) and satisfies (fl) and (/2) with a 
Lipschitz constant L < 4d. Then the system I.V.5J) /ias a unique solution. 

PROOF. We apply the classical result quoted in (E) of Section 12.21 to X = 
X* = R( n_1 ) , endowed with the Hilbert-Schmidt norm, Bu = fu — A £ u, 
A £ = XJ l D £ XJ f(u) = {f{ui)} and b £ X with components bi = -g n {^) - 

n d ¥{Di). The set of vectors of X, {UJ Up v i = l,--- , {n-l) d }, (Lfo)i = 

(pp. (j^J ) where i = 1, . . . , (n— l) d and j is the j-th vector of I d ordered with 

the lexicographic order, is an orthonormal system of X of eigenvectors of A £ 
with eigenvalues 

A/3 ^(Plc, 1+ ... + pjc, d ) Sin2 (^ 



a = - \ 2 ' , (3.9) 



satisfying 4f < c t < 1 for every 1 < I < n - 1. Thus, 4|/3| 2 < (A^ < vr 2 |/?| 2 . 
This property together with the obvious bound |\fr(e/3)| < ||*||£i(Rd) = 1, 
imply -A/j/%/3) > 4|/3| 2 > id. 

The operator B is strictly monotone. Indeed, for any u, v £ X, u ^ v, 
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we write u-v = E/3e/^( u /3 ~ v p){%) u p- Then, 
(Bu — Bv, u — v) 

= (~ T, TT^tW ~ "ft) U U a + (f(u)-f(v)),u-v) 



2 



+ ^ (Uj - Vi)(f(Ui) - f(Vi)) 

> Ad\u - v\ 2 - L\u - v\ 2 = (Ad - L)\u - v\ 2 , 

where we have applied ()2.3|) . 

Let us now prove that B is coercive. Since B is strictly monotone, and 
B(0) = /(0), we have 

(5u,u) = (B(u) - B(0),«) + (5(0),?!) > (4d- L)\u\ 2 - \(B(0),u)\ 
> [Ad-L)\u\ 2 - |/(0)| |4 

Therefore, 

>(4d_L)H- |/(0)|, 

\u\ 

which implies lim^i^^ ^ B ^j^ = oo. 

Finally, since for any u, v £ X, the functions t i— > U t D e U(u + to) and 
£ i — ► f{u + tv) are continuous, B is an hemicontinuous operator. 

□ 

3.3 Properties of the regularized and truncated kernels 

This section is devoted to prove some integrability properties of the kernels 
G e D n and estimates of the discrepancy between G £ D and Go, G £ D n , respec- 
tively. 

Along the section, 9 is a fixed positive real number satisfying > 2d — A. 

Lemma 9 Let e := e(n) = n ( 2d -^- s )/ 9 . There exists a constant C(6) > 
such that 

sup sup || G$i{x,-) \\ L 2 {D) < C(6). (3.10) 

n>l ieD 

Proof. The system {vp(n n (y))} is orthogonal in L 2 (D). Moreover, 
su Pa;eD \ v p( x )\ < C an d | Ajsl > 4|/3| 2 . Thus, 

II^D,n^>-)llL2 (D) - C ^ A 2 ^ — e- iTpTe- 
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for any > 0. 

The sum in the right-hand-side of this expression is comparable with a 
Riemann sum for the integral of the function |x| _4 ~ 9 on a region away from 
the origin. Observe that yd < \/3\ < nyd, for (3 S I d , and | 4 61 is radial 
and decreasing. Denoting by Qp = nf=i(/3« ~~ 1> A)> we have \[3\ > \x\, and 



1 



71" 



E 

pel* 



1 



4+6» 



< 2° 



< C 



(1/2,1)"* 



Thus, 



\\^D,n{ x '')\\L 2 (D) < —^~ n n - — 
Hence, choosing e = n^ 2d ~ 4: ~ e ^ e we obtain (|3.1()[) . 



-n 



2<2-4-0 



□ 



The next result gives an estimate in the L 2 (D x D)-norm of the approx- 
imation of the smoothed Green function by truncation. First, we recall the 
following facts (see Lemma 3.2 ^^): For any G I d , i,z£ 



\v p (x) - v p (z)\ <C\p\\x-z\ 



n*\/3\' 



1 



< 



C 



(3.11) 
(3.12) 



Lemma 10 Set e(n) = n^ 2d 4 9 ^ e . For every 7 > there exists a constant 
(7(7) such that 



sup\\Gp(x,-)-G^:(x,-)\\ LHD) < C{i)n 
xeD 



-,e(n) . 



(3.13) 



Proof. For simplicity, we write e instead of e(n). By the definitions of the 
kernels G e D and G £ D n , we can write 

/ \G £ D (x,y) -G s D<n (x,y)\ 2 dy < 2{A{x) + B(x)), 
Jd 

2 d V(e(3) 



with 



A{x) 
B(x) 



D 



E 

Pel*\I* 



-Vj3{x)vp(y) dy 



D 



E 

pel* 
2 d ^{e(3) 



-2 d ty(sf3) 



71 



Vf){x)vp{y) 



A 



P 



Vp(K n (x))v/3(K n (y)) 



dy 
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Since ^{e0) is a rapidly decreasing function, owing to the orthogonality of 
the system defined by the functions vg, we have 



£U ^ IAI4+6 



1 „ C{b)n° 



r d 1 dr 



4+6 



C(6)n 



(2d-4)(l-£) 



for every b > d — 4, where is a positive constant depending on b. 

Fix 7 > and choose b>2d-A> d- 4 such that (2d-4)(l-|) = -27. 
We obtain sup^g^ ^4(x) < C{^)n~ 2 "< . 

Clearly, < C{B 1 (x) + B 2 (x)), with 



Bi{x) 
B 2 {x) 



D 



D 



E 

Pel* 



m 2 



(v 3 (x)v a (y) - v 8 (n n (x))vg(K n (y))) 



dy, 



pel 

By virtue of (jTTTj) 



i>{£l3)v B {K n (x))v B {K n (y))) 



Bj{ x ) <c t ^ {£(3)l 1/31 < r nV * rd ~ ldr 



pel* 



\f3\ 2 n ne l A/2 



-i+z 



2d-l-2 



dy. 



for any / > 0. 

1 

Let I be such that 2d — 2 — (2(f — 4)^ = —7; we obtain sup zg£ , (x) < 
C(7)n-T. 

The orthogonality of the vectors ug(K n (x)), (j<3.12j) and the properties of 
^ imply 



B 2 {x) < C ]T 
pel* 



1 1 

+ 



A/3 ' vr2|/3|2 



i^)i 2 ^E|a_L_ 



< 



C{k)n d f n ^ r d - 1 dr C{k)n d d _ k _ 2 



n 2 e k 



1/2 



„2+fc 



< 



n 2 e k 



-n 



for any k > 0. Choosing A; such that (2d — 4)(l — |) = —27, we obtain 
su P:!;eDJ B 2 (x) <C( 7 )n- 2 ^. 

We have thus finished the proof. 



□ 
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Lemma 11 With the same assumptions as in Lemma^ set e(n) = 
n (2d-4-6)/9 _ There exists a constant C := C(X,6,d) not depending on n, 
such that 

\\Gd - G £ d ,Jl«(DxD) < Cn~\ (3.14) 
with 7 = |(0 + 4 -2d). 

Proof. It is a consequence of Lemmas 171 and UTil □ 
Remark In Lemma 111! the parameter can be chosen arbitrarily large. 
Therefore, for any given 5 > 0, one could obtain ()3.14)) with 7 = A — 5. On 
the other hand, we do not have an explicit control on the dependence of 
the constant C on 0. Actually, this constant appears in the formulation of 
the rapidly decreasing property of the Fourier transform of the regularising 
kernel ^f. 



3.4 Properties of the approximations 



One of the consequences of the properties of the modified Green kernels 
established in the preceding section is the following a priori estimate for the 
solution of equation 1)3.7)1 . This is the main result. 

Proposition 2 Fix > 2d — 4. Assume that the Lipschitz constant in 
(/2) satisfies L < min(4d, [C(0)]~ 1 ), with C{9) given in and the 

hypotheses of Proposition^ Let e(n) = n^ 2d ~ 4: ~ e ^ e \ Then, for any p G 
[I, a'] 



sup 

n> 



for some positive constant C{9) depending on 9. 

Proof. Let us write e instead of e(n). Since a < 2, Holder's inequality and 
()3.1U)) yield sup n sup x6 £, \\G £ D n \\L a (D) ^ C- Then, Holder's inequality and 
the properties on / imply 



K(*)l < 



D 



G £ D Jx,y)dF(y) 



+ SUp SUp || <?£),»» || 2,«(D) 



X (sup\\gn\\ La ' (D) +M + /(0) +L\\Un\\ La >( D )J ■ 

Notice that since u £ n is a step function, its L a -norm is finite. 

Proceeding as in ()2.5)) with the Green function Go replaced by G £ D n and 
using (j3.10)) . we see that for some positive constant C(9), 



E 



D 



G £ D J;y)dF(y) 



L a '(D) / 



< C(9). 



Since L < [C(0)] , the announced result follows easily. 



□ 
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3.5 Convergence results 

We devote this section to the proof of the approximation of the solution of 
(|1.1|) by means of the discretized scheme defined in (j3,7j) , with an appropriate 
choice of the smoothing parameter e. Here is the statement. 

Theorem 4 Fix A £ (0, 1) and 9 > 2d — 4. We assume that the hypotheses 
of Theorem and Propositions Q] and [H are satisfied. That is, if £ C a 



With a G (°» (d-2)(2-A)V(d-l)A / 

7j(6 + 4 — 2d). Then, there exist a constant C{^) > such that for any 



pe [l,a'], 



where e(n) = n 



A n ) I 



andL < min(a,4d,Cf 1 ,[C(^)]" 1 ). Set 7 

>r ai 

(3.15) 



" ~ u n"'\\Lv(n;L?(D)) < C(i)n 7 2^, 

e ), n> 1. 



Proof. Throughout this proof we write e instead of e(rt), for the sake of 
simplicity. Set 



T(x) 



+ 



D 



D 



D 



G D (x,y)-G £ Dyn (x,y)]f(u £ n (y))dy 

GD(x,y)[g(y) - g n (y)]dy+ jJG D (x,y) - G e D n (x,y)]g n (y) dy 

(3.16) 



D 



+ / [G D (x,y)-G%Jx,y)]dF(y), 



so that 

u{x) - u £ n (x) 



D 



G D (x, y)[f{u{y)) - f(u £ n (y))] dy + T(x). (3.17) 



We multiply both sides of (|3.17j) by f(u(x)) — f{u e n (x)), then we apply the 
inequality (|2.3j) and integrate over D and apply J222). We obt am 

—L\\u — u^||i2( D ) 



< 



+ 



( / G D (x, y)[f{u{y)) - /«(y))] dy) (f(u(x)) - /«(x))) dx 



id 

< -a 

+ 



D J D 

T(x)(f(u(x))-f(uUx)))dx 



D 



D 



G D (x,y)[f(u(y))-f(u £ n (y))] dy 
T(x)(f(u(x))-f(uUx)))dx. 



2 

L*(D) 



Now, by (1371711 . 



G D (x,2/)[/(«(y))-/«(y))]dy 



l n ~~ U nllz, 2 (,D) + 



I?{D) 

u(x) — u £ n (x))T(x) dx. 



D 
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Thus, 



-L\\u — it! 



e i|2 

n\\L 2 (D) 



< —a\\u — u' 



e 1 1 2 

n\\L?(p) 



a\\T\ 



L a (D) 



+2a (u(x)-u £ n (x))T(x)dx + (f (u(x)) - f (u £ n {x)))T (x) dx . 



D 



D 



From here, applying Holder inequality and the hypothesis (/l) and (/2) 
(a - L)\\u - u e n \\ 2 L 2 {D) 



<2a / (u(x)-u £ n (x))T(x)dx + / (/(«(*)) - f{u £ n (x)))T(x) dx 

JD JD 

< 2a\\u - u e n \\ La > {D) \\T\\ La{D) + ||/(u) - /«)|| iQ / (D) ||T|| La(D) 

< 2a\\u - u e n \\ La ^ D) \\T\\ La{D) + \\T\\ La{D) (2M + L\\u - u £ n \\ La ^ D) ) 

= (2a + L)\\u - u e n \\ La , {D) \\T\\ La{D) + 2M||T|| r « (£)) . (3.18) 

Notice that, by virtue of Lemma H3 and Proposition [2J for the choice 
e := e(n) and for any 1 < q < a', 

suvE{\\u-ul\\ q La , {D) )<C. (3.19) 

Assume p G [1,2]. Applying Holder's inequality, then (|3.18|) . Schwarz in- 
equality and (|3,19|) with q = 2 yields 

1/p , \ 1/4 



(e{\\u-uI\\1 2{d ^ < C(E(\\T\\1 HD) )} 



(3.20) 



Assume now p g]2, a']. Owing to ()3.18j) . Schwarz inequality and (|3.19|) 



E(h-<\\ p LHD) : 



2/p 



< C[E[ \\T\\ p La{D) 



(3.21) 



Hence, we have reduced the proof of the theorem to that of giving estimates 
of the L p (f2)-norm of H^ll^am), for p 6 [2, a']. 

By the definition of T(x) (see IjH.lfiJl ) and Holder's inequality, for any 
a > 1 we have 

H^Hl^d) < \\Gd — G £ Dtn \\ La ( Dx 0)11/(^)11 L a '{D) 

+ \\GD\\L a (DxD) \\9 — 9n\\ L c,>( D -) + \\G D ~ G e Dn \\ L c>( DxD ) 1 1 <?n 1 1 £a' (£)) 



+ 



D 



D 



dx 



1/a 



[G D (x,y)-Gh, n (x,y)]dF(y) 

Lemmas lll| |2] and the hypotheses on the coefficients / and g yield 

[|r|| L a (D) < C(7)n -7 

x [M+l+L 



\ 1/a \ 
\n £ n (y)\ a ' dyj + |/ 2 (0)| + \\g\\ La > iD) ) 

+ Cm" 1 + || J D (G D (;y) - G £ D>n (;y)dF(y)\\ La , {D) , (3.22) 
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where 7 = f (0 + 4 -2d). 

Let p G [1, a']. According to Q2.5JI . where we replace Gd by Gd — G^, n , 
we obtain 



J^(Gi?(-,y)-G!,, B (-,»)dF , (»)||J a , (J3) )V P 
<c(J d \\G D (x, •) - G s D , n (x, Wl^d) dx 



l/a' 



= C^J d \\G D (x, •) - Gh, n (x, -m^WGDix, ■) - Gh >n (x, 
<C\\G D (x r )-G% )n (x,-)\\%^ DxD) , 



l/a' 

\l<*(d) d x 



by Lemma |21 and (|3.10j) . Together with (|H.22jl . Proposition |21 and Lemma 
fTTl this implies 

(EiWTWUD))) 1 P < C(7) (n^ + n-^) < C( 7 )n^. 

Substituting this bound in IJM.20J) . l|3.21j) yields the upper estimate (j3.15jl 
and therefore, the theorem is proved. 

□ 

Remark: In the previous theorem, the values of p cannot be arbitrarily 
large; therefore, one cannot obtain a rate of convergence for the sample 
paths, as is the case for d < 3 (see for instance Corollary 2.4 in 
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